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In this paper, we present measurements of the ortho-positronium emission energy in vacuum from 
mesoporous films using the time of flight technique. We show evidence of quantum mechanical 
confinement in the mesopores that defines the minimal energy of the emitted Ps. Two samples 
with different effective pore sizes, measured with positron annihilation lifetime spectroscopy, are 
compared for the data collected in the temperature range 50-400 K. The sample with smaller pore 
size exhibits a higher minimal energy (73±5 meV), compared to the sample with bigger pores (48±5 
meV), due to the stronger confinement. The dependence of the emission energy with the temperature 
of the target is modeled as ortho-positronium being confined in rectangular boxes in thermodynamic 
equilibrium with the sample. We also measured that the yield of positronium emitted in vacuum is 
not affected by the temperature of the target. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Positronium (Ps), the bound state of electron and 
positron, was extensively investigated since its discovery 
1] contributing to the development of bound state QED 
(see e.g. 0] for a review on the current status of this 
field). Furthermore, this system provided stringent lim- 
its on possible deviation from the Standard Model that 
could indicate new physics Q (see also [|| for a compre- 
hensive review of former experiments). Moreover, in the 
field of materials science Ps found various applications 
due to its unique properties (see e.g. [5, 6] for mod- 
ern reviews on this subject). One recent example is the 
characterization via positron annihilation lifetime spec- 
troscopy (PALS) of low-k dielectrics that are potential 
candidates for the next generation of integrated circuits 
01- 

The motivation of the work presented in this paper is 



to understand if mesoporous silica with an interconnected 
pore network could be used for producing a high fraction 
of positronium at low temperatures. This would open the 
door for a new generation of experiments in fundamen- 
tal research. Cold positronium could be used to improve 
the precision of spectroscopic studies of Ps [t| [l(| or to 
perform the first spectroscopy of the Ps 2 molecule fTTj j . 
Furthermore, it could provide an alternative to the meth- 
ods that are used for anti-hydro gen formation 12|-[l4j|. 
As it was suggested sometime ago 15|-[lj}, anti- hydrogen 
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could be formed using charge-exchange of Ps with anti- 
protons. This was demonstrated for the charge conjugate 
reaction [l8[ and in the ATRAP collaboration resonant 
charge exchange collisions of positrons with Rydberg Cs 
atoms were used to form Rydberg Ps that via charge ex- 
change with the anti-protons produced anti-hydrogen in 
Rydberg states [lj| . Recently, two experiments (20, Hl| 
were proposed to perform an anti-gravity test using this 
process. In both experiments, one of the main issues is 
that a high fraction of Ps at low temperatures should be 
available. Another interesting application would be the 
possibility to perform an experiment in order to confirm 
the interpretation of the recent DAM A/LIBRA annual 
modulation signal p2l [23j | as generated by Mirror type 
dark matter |24|. A step further would be to achieve 
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Bose-Einstein condensation of Ps [25[. This would al- 
low the exploration for the first time of the effects of the 
collective properties of a matter-antimatter system. 

In mesoporous silica, Ps is produced by injecting 
positrons into the film and the distribution of the implan- 
tation depth follows a Makhovian profile [26[ . In the fol- 
lowing we solely consider the long lived triplet spin state 
(called ortho-positronium with 142 ns lifetime, o-Ps) be- 
cause the singlet spin state (called para-positronium, p- 
Ps) has a very short lifetime of 125 ps and can be consid- 
ered as annihilating in the target. The o-Ps (for simplic- 
ity we will refer to it as Ps) that diffuses into the pores 
loses its kinetic energy via scattering. If the pores are in- 
terconnected, the Ps has a probability to tunnel from one 
pore to another. A fraction of the Ps reaches the film sur- 
face and exits into vacuum. A classical model of thermal- 
ization process was developed by Nagashima et al. [27j . 
Their calculations reproduce very well the behavior for 
SiC>2 aerogel with pore sizes of about 100 nm. However, a 
classical approach is not expected to give reliable predic- 
tions for Ps confined in few nm pores because quantum 
mechanical effects become relevant. As a matter of fact, 
in this regime the de Broglie thermal wavelength of Ps is 
comparable with the size of the pores. Recently, Mariazzi 
et al. [2S] considered phonon scattering to reproduce the 
thermalization process of Ps in a box (closed porosity) 
showing that the minimal energy is not that of the low- 
est accessible level because the momentum phonons can 
exchange is fixed. In the same paper, they pointed out 
that this is not the case in rectangular channels because 
in one direction the side of the potential well tends to in- 
finity (z-axis). Therefore, the magnitudes of the k x and 
k y momentum are quantized but the k z tends to be a 
continuum and the minimal energy that Ps can reach is 
given by the ground state in the x-y components. Thus 
measuring the Ps emission energy provides a method to 
distinguish between the two different pores architectures 
(see Section HVCl) . 

In previous studies of Ps emission in vacuum using 
time of flight (TOF), many interesting effects, such as 
the emission from the surface of different materials [29| - 
3l| were investigated. Recently, this technique was ap - 
plied to study mesoporous and hybrid silica films |33 |- 
[37j , to evaluate the continuity barrier [38| and the effect 
of thermalization for pore surfaces decorated with dif- 
ferent groups [lH . However, the influence of the tem- 
perature on Ps emission in vacuum from mesoporous 
thin films was never studied in detail. To our knowl- 
edge, only a work of Mills et al. [3l| with fumed silica 
revealed a dependence with the temperature of the Ps 
emission in vacuum. The fraction of Ps in the low en- 
ergy tail was estimated but no quantitative estimate on 
the value of the minimal emission energy was obtained. 
Very recently Cassidy et al. [Io| measured the emission 
energy of Ps from mesoporous films using Doppler spec- 
troscopy. This is a different technique from that used 
in the present study. Results consistent with the one 
presented in this paper (though not sample temperature 



dependencies) have been obtained. 



II. DESCRIPTION OF THE EXPERIMENTAL 
APPARATUS 



A. Positronium production 

In this paper, we study the Ps yield and emission 
energy in vacuum as a function of the target temper- 
ature for two different kinds of mesoporous thin films 
with the same tetraethoxysilane (TEOS) mineral source 
for the silica network skeleton precursor: CTAC1-TEOS 
and F127-TEOS. The density of the C sample is approxi- 
mately 1.2 g/cm 3 and of the F sample is 1.5 g/cm 3 . Both 
samples were spin-coated on glass similar to the ones 
we measured in 



411 ] . The C samples are prepared via a 



sol-gel process using cetyl trimethyl ammonium chloride 
(CTAC1) cationic surfactants as the organic pore gener- 
ator (porogen) agent [42j]. A pure aqueous method is 
used. The CTAC1/TEOS molar ratio for the films pre- 
pared is 0.22. After deposition, the CTACl-TEOS/Glass 
samples are treated at 130 °C and stored in air. The F 
samples use non-ionic Pluronic F-127 triblock copolymer 
(EO106PO70EO106) as surfactant and were prepared in 
the same way as described in 43(. Both samples were 



calcinated for 15 minutes at 450 °C in air immediately 
before the e + measurements. The recorded X-ray diffrac- 
tion patterns indicate no symmetry in the pore organi- 
zation. 



B. Slow positron beam 

The ETHZ slow positron beam used for these measure- 
ments is described in greater detail in [44| . The positrons 
flux is 25000 e + /s. The slow positron beam is stopped in 
the SiC>2 target. The positrons can either form positro- 
nium, i.e. o-Ps or p-Ps, or annihilate into 27's. The 
detection with an MCP (Micro-channel Plate) of the sec- 
ondary electrons (SE) emitted when the positrons hit the 
target serves for tagging the positronium formation. The 
SE leave the target accelerated to 1-11 keV by the same 
voltage applied to the target relative to the grounded 
transport tube that is used to implant the positrons in 
the positronium converter. The SE are then transported 
by a magnetic field in the backward direction, as shown 
in Figure [T] The electrons move along the magnetic field 
line in spirals and are deflected to the MCP region by 
the E x B filter. The tagging efficiency varies from 70% 
to 30% in the energy range 1-10 keV. 

The samples are mounted on a cryocooler head to allow 
the possibility of varying their temperature in the range 
of 50-400 K. 
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FIG. 1: (Color online) Experimental setup for the TOF mea- 
surements (the scale is in mm). The dashed line is the tra- 
jectory of the incoming positrons (blue) and the dotted line is 
the one of the secondary electrons (red). 



C. PALS and TOF detectors 



III. MONTE CARLO SIMULATIONS 

In order to design the detectors and interpret the data, 
simulations served as a powerful tool. In the Monte Carlo 
(MC) simulations of our setup, the 3D EB-fields were cal- 
culated with the COMSOL multi-physics program and 
the positron/electron trajectories in the beam were sim- 
ulated with GEANT4. The simulation of the photon de- 
tection in the apparatus was based on the same package 
(47} . New classes were written in order to simulate the o- 
Ps production, propagation in the beam pipe, reflection 
on pipe walls that were assumed to be Knudsen-like (Ps 
is reflected isotropically) , pick-off effect or decay (Ps is 
not a standard particle in GEANT4). The events for the 
o-Ps— > 37 process were generated taking into account the 
decay matrix element [481 ] . The geometries of the beam 
transport pipe, photon detector, positron tagging sys- 
tem and its material were coded into simulations. The 
results were cross-checked with our experimental mea- 
surements for both photon detection [J [4(| H^, [5(| and 
particle transport in the EB-fields [44l ]. 



The start time to for the detectors is triggered by the 
MCP detecting the SE emitted when the positrons hit 
the target. In both PALS and TOF detectors the stop 
is given by one (or more) annihilation photons deposit- 
ing some energy in the calorimeter (ECAL). Both ECAL 
are composed of BGO crystals with hexagonal shape, 61 
mm external diameter and 200 mm length. The time res- 
olution of the system MCP-ECAL was measured to be 
around 5 ns (FWHM). The typical energy resolution of 
the crystals is about 25-30% (FWHM). The ECAL for 
the TOF is placed behind a lead slit at a distance z from 
the target that can be varied (as sketched in Figure [TJ. 
The detector is screened by four half cylinders of lead 
surrounding the beam pipe. The thickness of the shield- 
ing is 70 mm, the width of the slit is set to 5 mm and 
its position with respect to the target can be adjusted. 
For this measurement in order to maximize the signal to 
background ratio the center of the slit was placed at a 
distance of 18 mm from the target. The main contribu- 
tion to the background is given by photons coming from 
direct positrons and p-Ps annihilations in the target (so 
called prompt peak). Some of those photons can be de- 
tected in the crystals after Compton scattering in the 
lead shield. The determination of the slit position with 
respect to the target was done by scanning the slit posi- 
tion in 0.1 mm steps and recording the maximum of the 
511 keV annihilation peak. The scans were performed 30 
min after the temperature was set on the sample. This 
technique was very important in order to correct the slit 
position as a function of the temperature. In fact, at 50 
K the contraction of the cryocooler head was measured 
to be 1.30 ± 0.05 mm in agreement with the prediction 
of finite elements calculations performed with the COM- 
SOL package [45) . The PALS detector was designed to 
have a large acceptance to provide a uniform efficiency 
for detecting the Ps emitted in vacuum |46| . 



IV. RESULTS 

The measurements were taken in a clean vacuum of 
10 -9 mbar. To avoid water contamination of the film 
during the cooling down, the target was kept at room 
temperature for an hour using a heater before lower- 
ing its temperature. The cooling cycles were repeated 
and the data confirmed the reproducibility of the results. 
The signals from the BGO's photomultipliers are split to 
record both energy and timing with a charge-to-digital 
converter (QDC CAEN v792) and a time-to-digital con- 
verter (TDC CAEN v775). A cut on the energy deposited 
in the BGO between 300 < Ebgo < 550 keV was applied 
to optimize the signal to background ratio suppressing 
Compton scattering events in the collimator from direct 
and p-Ps annihilations in the target. 

A. PALS measurements 

In this section we present the results we obtained from 
the PALS measurements. The spectra are analyzed us- 
ing the LT9 [52j program. The lifetimes of the decay 
components and its fractions are resolved by fitting the 
PALS spectra. The program finds 3 exponentials con- 
voluted with a 5 ns FWHM resolution function of the 
spectrometer. The shortest exponential (less than 4 ns) 
is originated by direct positron and p-Ps annihilations. 
It is disregarded because we are interested only in con- 
tribution of o-Ps. We define as (r 2 , h) ,and (t v , I 3 ) the 
intensities and the lifetimes of the two longer exponen- 
tials. To determine the yield Y v of Ps emitted in vacuum 
and the lifetime Tf in the pores of the film we used a 
model of Ps escaping in vacuum 51]. According to this 
model the lifetime ry in the pores of the film is defined 
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Tf = ({T 2 - 1 -T v - 1 )h/{h+h)+T v - 1 ) (I) 

The yield Y v of Ps emitted in vacuum was calculated 
according to: 



Y v = {h + iaK/Cr/- 1 + k v ) 



where 



(r 2 - 



^)h/{h+h) 



(2) 



(3) 



is the escape rate of free Ps into vacuum. We extracted 
the films thicknesses for the two samples by fitting the 
total Ps yield as a function of the positron implantation 
energy to the Makhovian profile. These thicknesses are 
reported in Table U in which we also present the values 
of Tf and k v calculated with the expressions above. We 
use the results of the fits of the PALS spectra at 6 keV 
for the C sample and fO keV for the F sample at 50 
K and 300 K. These implantation energies will serve as 
a reference for the rest of the paper. We choose those 
values in order to maximize the amount of thermalized 
Ps. At these energies the majority of the positrons are 
still implanted within the films (no significant drop of the 
total Ps yield I tot is observed, see Fig. [5]) and the emis- 
sion energy is in the constant region (see next Section). 
With Tf one can calculate the effective pore size a in the 
films applying the Gidley et al. @, [53| extension of the 
Tau-Eldrup model [EIj]-[57| (we will call it hereafter RTE 
model). As one can see, at lower temperature the life- 
time in the film increases as predicted by the RTE model. 
This can be understood in the following way: the over- 
lap of the Ps wave function with the volume contained 
within a distance 5 = 0.18 nm 0, [53| from the walls for 
which the annihilation rate is assumed to increase is less 
for Ps confined in the pores occupying the ground state 
than for Ps in excited states. Since at lower tempera- 
tures the population of the ground state is higher the 
pick-off rate decreases. However, the measured lifetimes 
are lower than what is expected by the calculations us- 
ing the RTE model. For the C sample the discrepancy is 
less than f0% but for the F sample the measured value 
is 30% lower suggesting that the pick-off rate increases. 
Deviations from the RTE model were already observed in 
previous measurements using the 2 to 37 ratio technique 
|58l [50( 1 and in a measurement using PALS [6(| ■ Differ- 
ent factors could be responsible for that as pointed out 
in those papers and deserves further studies. 

Since the parameter S of the RTE model was calibrated 
at room temperature we use the lifetimes at 300 K to de- 
termine the effective pore size that Ps experiences. We 
report the results at the end of the paper in Table |TT] 
for both rectangular channels and cubic boxes. Interest- 
ingly, we found that the yield of Ps emitted in vacuum 
can be considered as independent of the temperature of 
the target (see Fig. [2]). This may seem to contradict 
the fact that the lifetime in the films increases with the 
temperature. However, our measurements indicate that 
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FIG. 2: Yield of Ps emitted in vacuum (Y v ) and total yield of 
Ps (hot = I2 + h ) for F (upper plot) and C samples (lower- 
plot) at room temperature and 50 K. 



TABLE I: Film thickness (Z), Ps lifetime in the pores Tf and 
escape rate k v at 50 K and 300 K. 



Sample 


Z[nm] 


T f 0K [ns] 


T^ ox [ns] 






C 


700±200 


54±1 


60±1 


27±1 


25±1 


F 


1000±200 


74±1 


82±1 


17±1 


13±1 



the escape rate in vacuum is smaller at lower tempera- 
ture (see Table |T| , explaining the observation that the 
yield is constant. A possible explanation of this effect 
can be found considering that Ps tunnels from one pore 
to the other. In this case, the tunneling probability will 
decrease with temperature explaining our measurements. 



B. TOF measurements 

Figure [3] shows the acquired TOF spectra at different 
implantation energies for the F sample. For every run 
the time spectra of each crystal are calibrated finding 
the position of the peak arising from the annihilations 
in the target that defines to. First we extract the mean 
energy of Ps emission in vacuum applying the analysis 
method proposed in [3l| as follows. The background due 
to the annihilations in the target is subtracted by fit- 
ting the measured time spectra with the resolution func- 
tion of our detector determined using targets (aluminum 
and kapton) in which the Ps formation is negligible. The 
TOF spectra are corrected for the Ps decays and the time 
spent in front of the detector with the factor l/t- e +t / T ". 
The maximum of the peak distribution defines the mean 
Ps emission energy in the direction perpendicular to the 
film surface. Let us note that, neglecting the reflection of 
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FIG. 3: TOF spectra of the F sample for positron implantation 
energies of 0.7-1-4-10 keV. 



Ps in the beam pipe, with the time-of-flight method one 
measures only the mean of the energy component per- 
pendicular to the collimator that we define as the z-axis 
(we define it as < E z >). The triangles in Fig. HI rep- 
resent it as a function of the implantation energy. One 
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can see that starting from 3 keV for the C sample (this 
was confirmed by other measurements on similar sam- 
ples [13]) and 4 keV for the F sample the value of the 
energy emitted in vacuum tends to be constant. Clearly, 
the emission energy of Ps calculated in this way is not 
the minimal energy due to the confinement in the pores 
because one has different contributions given by the con- 
volution of the emission energy with the implantation 
profile. In order to isolate the thermalized part from 
the TOF spectra, one has to subtract the contributions 
of non-thermalized Ps components (see Fig. [5]). To es- 
timate these non-thermalized contributions, we suppose 
that the shape of their distribution, NT(t), is represented 
by the TOF spectrum obtained at a low implantation en- 
ergy. To select the implantation energy (from the ones 
we had measured) and the scaling factor of NT(t) for the 
subtraction, we relied on the MC. We used the values for 
which the best fit between the MC and the spectra ob- 
tained after the subtraction of NT(t) was achieved (for 
more details see [HI]). We found that the best fits were 
obtained for 2 keV with L2fceV=200 nm in the case of the 
C sample, and 3 keV and L 3 £, e ^=350 nm for the F sam- 
ple. For a given implantation energy E,, NT(t) is scaled 
down by the fraction of positrons implanted at depths 
smaller than L 2 fc e y and L 3 fc e y. This was determined by 
using a Makhovian profile. 



data: CTACI-TEOS 300K-6kV 

o Target component 

+ Target + 2kV correction 
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FIG. 5: Upper plot: the solid line is the TOF spectra without 
correction, the diamonds is the contribution from the target 
and the crosses is the sum of the contribution from the target 
and the non-thermalized part. Lower plot: the TOF spectra 
after the correction of the target and the non-thermalized com- 
ponents. 



FIG. 4: Positronium mean emission energy < E z > as a 
function of the implantation voltage at a target temperature 
of 300 K. The triangles represents the energy extracted with 
the maximum analysis method (MA), the squares are after the 
subtraction of the non-thermalized part (maximum analysis 
method after correction of the spectra, MA-CS). 



The results of this analysis are shown as the squares 
in Fig. 01 A parabolic fit is used to determine the posi- 
tion of the maximum. The statistical error of the fit is 
typically ±9 ns for the F and ±6 ns for the C sample. 
The uncertainty on the determination of the slit position 
of ±0.1 mm results in a systematic error of the order of 



6 



±1 meV in the determination of the Ps mean emission 
energy < E z > (for implantation energies above 3 keV). 
The subtraction procedure described above introduces a 
systematic error that we estimated analyzing the data 
using different values of L*2keV ± 50 nm and h^keV ± 50 
ran. The estimated error is ±2.2 meV for the C and ±1.6 
meV for the F sample. Thus the combined statistical and 
systematic error is at a level of ±2.9 meV for the F and 
±3.0 meV for the C sample. 

The results for the C and the F samples at room tem- 
perature and at 50 K are shown in Figs. [S][7J For implan- 
tation energies higher than 4 keV for the C and 5 keV 
for the F sample the values < E z > of the mean emission 
energy are constant. In the C sample, < E z > reaches 
its constant value at lower implantation voltages because 
the pore size is smaller than in the F sample (see Fig. [B]). 
Those values are higher than the thermal energy that Ps 
will have if it would thermalize at the temperature of 
the film. As expected in the presence of confinement in 
the pores, the mean emission energy is higher for the C 
sample with pore sizes smaller than the F sample. 
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FIG. 7: Upper plot: Ps mean energy < E z > for implantation 
energies higher than 2 keV at 50 K and 300 K for the C 
sample. Lower plot: Ps mean energy as a function of the 
implantation energies higher than 4 keV at 50 K and 300 K 
for the F sample. 
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FIG. 6: Ps mean emission energy < E z > as a function of 
the positron implantation energy for C and F samples at 300 
K. 



The TOF technique measures < E z >, the mean en- 
ergy of the Ps atoms in the z-direction. To find the mean 
emission energy of Ps in vacuum, < E z > should be mul- 
tiplied by a factor £ that takes into account the angular 
distribution. Assuming that the Ps is emitted mono- 
energetically and isotropically from the surface, one can 
calculate that £ is equal 2. In this estimation, the reflec- 
tion of Ps in the beam pipe and the detector acceptance, 
i.e. the fact that a fraction of events decaying before or 
after the collimator aperture are detected, are not taken 
into account. Therefore, to determine £ considering these 
effects we used the MC simulation we described in Sec- 
tion IIIII As shown in the upper plot of Fig. [5J a satis- 
factory agreement between the data and the MC (adding 
the spectra at 2 keV that takes into account the non- 
thermalized Ps) is achieved. We attribute the difference 
between 40 and 100 ns to the approximation used in the 
subtraction method where the contribution of the non- 



thermalized Ps is underestimated since only a spectrum 
of a defined energy is used for this correction. The fact 
that the Ps is emitted with an angular spread is clearly 
supported by the data. As one can see in the lower plot 
of Fig. [8] for Ps emitted with no angular spread the data 
are not reproduced. The physical interpretation is that in 
the films we studied the pores have no organization thus 
they are expected to be randomly aligned. The value of 
£ estimated with the MC is 1.7. This is consistent with 
the expectation of the analytical result and the values 
reported in previous experiments [3(J HH, [H, [H, Ho| . 

A detailed scan shows that the mean Ps energy (< 
Eps >= £ < E z >) decreases with the sample temper- 
ature down to a minimum level (see Fig. [9]). For the 
C sample this value is basically constant (73 ±5 meV) 
in the range of temperature in which we performed our 
measurement. This can be understood by the fact that in 
this sample the confinement energy is much higher than 
the thermal energy at room temperature (kT ~ 25 meV) 
thus almost all Ps is in the ground state. For the F sam- 
ple there is a weak dependence on the temperature. Due 
to the bigger pore size compared to the C sample, the en- 
ergy of the ground state is only twice the thermal energy 
at room temperature. Therefore, the probability to find 
the Ps occupying an excited state is higher. As expected, 
this probability decreases with the temperature thus the 
minimal energy reaches its constant value of 48±5 meV. 

The time that Ps spends in the films before being emit- 
ted in vacuum was not considered in our determination 
of the emission energy. The measurements presented in 
Fig. [5] are in a regime in which a classical approach is 
not expected to give reliable results. Some theoretical 
work to develop a full quantum mechanical picture of 
the emission process is required to address this problem 
(as pointed out in as well). 
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data: CTACI-TEOS 300K-6kV 



Monte-Carlo + 2kV correction: 
<E o Ps >=75 [meV] 





data: CTACI-TEOS 300K-6kV 
Monte-Carlo + 2kV correction 
E .p a =75 [meV]- no angular spread 



FIG. 8: Upper plot: comparison between the data of the C 
sample at 6 keV and the MC simulating mono-energetic Ps 
emitted isotropically from the film surface. Lower plot: com- 
parison between the data of the C sample at 6 keV and the MC 
simulating mono-energetic Ps emitted perpendicular from the 
film surface. In both cases, the measured non thermalized part 
(called 2 keV correction in the legend) was added to the MC. 



C. Discussion 

To understand the behavior of the value of the minimal 
energy as a function of the film temperature we present 
a simple model of Ps in thermodynamic equilibrium at a 
temperature T in rectangular boxes. 

The expectation value < H > of the Hamiltonian op- 
erator for Ps confined in 1 dimensional infinite well in 
contact with a reservoir at a temperature T is given by: 



< H >= kT- 



1 dZ 
ZdT 



where Z is the partition function defined as 



Z(a) 



(4) 



(5) 



where a is the dimension of the well, m is the Ps mass, 
n is the principal quantum number and h and k are the 
Planck and the Boltzmann constants. To calculate the 
mean value < H > of the energy for the 3D case we can 
use 



< H >=< H x > + < H v > + < H z > 



(6) 



where to calculate < H y > and < H z > one can substi- 
tute the pore side length a in Eq. [3] with b and c. For Ps 
confined in rectangular pores, we thus obtain: 



< H >= kT 1 



1 dZ{a) 1 dZ(b) 1 dZ(c) 



Z(a) dT Z(b) dT Z(c) dT 



■ F127-TEOS 

▼ CTACI-TEOS 

RECTANGULAR PORES 

CUBIC BOX PORES 




(7) 



100 150 200 250 300 350 

Sample Temperature [K] 



FIG. 9: (Color online) Positronium mean energy as a func- 
tion of the mesoporous film temperature. Those results are 
obtained at 6 keV for the C and 10 keV for the F sample. 
The solid lines are the results of the using Eq. Q] with the pore 
side lengths a,b,c left as free parameters. The dashed lines 
were obtained fitting with Eq. with a single side length free 
a — b = c (cubic box pores). 



For the case of a cubic box one can set < H x >=< 
H y >=< H z >. 

To compare the prediction of the pore size that one 
can extract from the TOF measurements with the PALS 
results we fit the data using Eq. [7] To construct the 
function used for the fit we kept only the first 50 terms 
of the sum. This is very conservative, the probability 
for Ps to occupy a state higher than n>10 for the kind 
of target and the temperatures we used in this study is 
already negligible. The solid lines in Fig. |H] represent 
the fit to the data where the pore side length (a,b,c) are 
left as free parameters. The fits were repeated assuming 
cubic box pores and the results are shown as the dashed 
lines in Fig_^ [3 We used MIGRAD from the MINUIT 
package [62[ as a minimization procedure. As one can 
see, the fit to the data suggests that the pores of both 
samples are better modeled as rectangular pores than 
cubic boxes. As proposed in [iol |. to compare the values 
obtained from the fit with the ones extracted from the 
PALS measurements, one has to add twice the parameter 
S (see previous Section). The pore side lengths obtained 
in this way are reported in Table ITT] In particular, the 
fit supports the idea that the pores are better modeled 
by rectangular channels. The pore side length in one 
direction (c) obtained from the fit is much longer than 
the substrate thickness. Since, the same values for the 
side length a and b are obtained fitting the data assuming 
rectangular channels instead of rectangular boxes we do 
not report this value in Table UU as it has no physical 
meaning. 

According to the quantum mechanical model for Ps 
thermalization of Mariazzi ct al. [28] the index n in the 
sum of Eq. [5j could differ from 1 for cubic pores if the 
level separation of two close Ps energy levels is higher 
than the maximum momentum that a single phonon can 
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TABLE II: Comparison of the pore sizes obtained from the 
PALS measurements at 300 K and from the fit of the Ps 
mean emission as a function of the temperature of the sample 
(see Fig. O for both cubic box (BOX) and rectangular pores 
(RECT). Minimal energy of Ps < Eps > for which the errors 
are the combined statistical and systematic error. 

C Sample 

TOF PALS 

o-box [urn] 3.3±0.1 4.2±0.5 

(a, 6) rect [nm] (2.7±0.8, 2.7±0.8) (3.3±0.5, 3.3±0.5) 

< E Ps > [meV] 73±5 

F Sample 

TOF PALS 

a B ox [nm] 4.1±0.1 6.4±0.5 

(a, 6) rect [nm] (2.9±0.8, 3.6±1.8) (4.3±0.5, 5.3±0.5) 

< E Ps > [meV] 48±5 



exchange. We performed fits with different n > 1 but 
the results did not improve, supporting the idea that the 
pores are not very well modeled by cubic boxes. To sum- 
marize, for the C sample with the TOF data the best 
fit was obtained for pores modeled as square channels of 
2.7±0.8 nm side length while the lifetime method gives 
3.3±0.5 nm. The best fit to the data for the F sample was 
obtained for rectangular channels with a cross section of 
2.9±0.8 nm by 3.6±1.8 nm that has to be compared with 
4.3±0.5 by 5.3±0.5 of the PALS measurement. This pore 
size was extracted by applying the RTE model to repro- 
duce the measured lifetime of 74 ns assuming the ratio 
a/b to be the same as in the TOF results. Both values 
obtained with the TOF measurement are systematically 
lower than the PALS results. Nevertheless, considering 
the approximation used in our model, the assumptions 
that the pores can be treated as rectangular channels 
and the uncertainty in the determination of the pore size 
we conclude that our results (summarized in Table [TTJ) 
are in reasonable agreement with expectations. 



V. CONCLUSIONS 

In this paper, we show that the yield of Ps emitted 
in vacuum measured with the PALS technique is inde- 
pendent of the temperature of the mesoporous thin films. 
The lifetime in the films increases with a decrease in tem- 
perature but the measured values are lower than expected 
by the RTE model. This suggests that another source of 
pick-off which depends on the temperature should be in- 
voked. Further studies are needed to investigate the ori- 
gin of this effect. The escape rate in vacuum decreases 
explaining the observation that the yield of Ps emitted 
in vacuum is the same at 50 and 300 K. 

Furthermore, we show that due to quantum mechani- 
cal confinement in the pores the Ps emission energy into 
vacuum has a minimal value. The minimal energy is 
higher for the sample with smaller pores and this con- 
stant value is reached at a lower positron implantation 
energy. Our results arc in fair agreement with a model 
of Ps confined in rectangular channels in thermal equi- 
librium with the sample. The measured minimal energy 
for the C sample was found to be 73±5 meV while for 
the F sample it is 48±5 meV. These experimental results 
provide a solid ground to understand how to produce Ps 
at lower temperature and could serve to develop a more 
realistic model to interpret the data. 
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